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1.  INTRODUCTION 


Preliminary  sizing  analyses  of  aerospace  vehicles  which  utilize 
solid  propellant  nocket  motors  usually  require  rather  extensive  trade¬ 
off  analyses  to  ensure  the  near-optimal  selection  of  pertinent  motor 
design  parameters.  For  multistage  vehicles,  the  motor  design 
analyses  usually  follow  energy  management  studies  (e.  g.  ,  Ref.  1)  which 
apportion  propellant  masses  optimally  among  the  constituent  stages  to 
achieve  prescribed  vehicle  performance  characteristics;  e.  g.  ,  to 
minimize  gross  weight  for  a  required  velocity  increment  or  to  maximize 
velocity  for  a  fixed  gross  weight.  The  objective  of  this  investigation 
was  to  develop  a  digital  computer  program  which  would  permit  rapid, 
accurate,  preliminary  design  tradeoff  analyses  of  solid  propellant 
rocket  motor  concepts. 


2.  PROGRAM  CAPABILITIES 


The  program  described  in  this  report  computes  rocket  motor 
length,  propellant  and  inert  component  weights,  mass  fraction,  burn 
time,  specific  impulse,  stage  burnout  velocity,  and  other  performance 
parameters  as  a  function  of  input  motor  diameter,  chamber  pressure, 
thrust,  payload  mass,  propellant  ballistic  properties,  propellant  web 
fraction,  port-to-throat  area  ratio,  residual  sliver  fraction,  and  other 
required  parameters. 

The  generic  representation  of  a  selected  propellant  grain  con¬ 
figuration  by  its  web  fraction,  Fw.  and  residual  sliver  fraction,  Fp,  and 
a  prescribed  port-to-throat  ratio,  Apt,  permits  evaluation  of  the  applica¬ 
bility  of  many  candidate  grain  designs  without  requiring  detailed  surface - 
web  calculations.  This  feature  simplifies  input  preparation  and  enhances 
program  flexibility  for  its  intended  use  in  conceptual  design  tradeoff 
analyses.  'Realistic  values  of  Fw  and  Fp  may  be  derived  from  known 
geometric  characteristics  of  previously  configured  grain  cross  sections 
(Refs.  2  and  3).  All  star,  wagonwheel,  cylindrical  port,  and  slotted 
tube  grain  configurations  wherein  the  propellant  perforations  traverse 
the  entire  cylindrical  motor  length  are  ideally  suited  for  analysis  using 
this  technique.  Tapered  grain  configurations  and  those  which  employ 
radial  slots,  or  longitudinal  slots  which  traverse  only  a  fraction  of  the 
cylindrical'  motor  length,  may  be  analyzed  provided  a  length-averaged; 
propellant  web  fraction  is  entered,  Ellipsoidal  head-end  propellant 
webs  and  inert  propellant  slivers  may  be  included  in  the  design  analysis 
simply  through  the  proper  selection  of  input  code  words,  as  described 
later. 

Dimensions  and/or  weights  of  inert  components,  including  the 
motor  case,  head  closures,  nozzle,  attachment  skirts,  interstage 


structure,  liner,  insulation,  and  igniter,  are  computed  using  proven 
theoretical  and  empirical  relationships  in  conjunction  with  input 
material  physical  properties.  Either  conical  or  contoured  nozzle 
designs  may  be  analyzed.  The  nozzle  expansion  ratio  is  sized  to 
provide  as  near  optimum  expansion  as  possible  while  maintaining  the 
nozzle  exit  diameter  less  than  the  input  maximum  allowable  value. 
(Optimum  expansion  occurs  when  the  nozzle  exit  pressure  is  equivalent 
to  the  input  ambient  pressure.  }  The  validity  of  the  computer  program 
was  confirmed  through  the  analytical  reproduction  of  various  rocket 
motor  designs  which  are  documented  in  Reference  4. 


3.  DESIGN  ANALYSIS 

3.  1  INTERNAL  BALLISTICS 

The  internal  ballistics  analysis  described  in  this  section  is 
based  on  the  assumption  of  constant  mass  flow  rate  throughout  motor 
operation.  The  appropriate  gas  dynamic  relationships  for  combustion 
gas  flow  within  the  motor  port  and  nozzle  are  based  on  the  assumption 
of  one-dimensional  flow  of  a  perfect  gas.  Most  of  the  applicable 

thermodynamic  relationships  may  be  found  in  standard  propulsion 
textbooks;  for  example,  Reference  5.  Following  program  initialization 
and  the  reading  of  appropriate  input  parameters,  the  rocket  mqtor 
design  analysis  proceeds  as  described  in  the  following  paragraphs. 

The  required  thickness  of  the  motor  case  wall  is  sized  using  the 
well-known  hoop  stress  relationship 

Pc  Dm  Fsm  . 

Tm  - — - —  .  w-i/ 

2crm 

The  independent  variables  Pc,  Dm,  and  Fsrn  are  program  inputs.  The 
yield  strength  of  the  motor  case,  crm,  and  most  of  the  other  required 
inert  component  physical' properties  discussed  later  are  specified 
within  the  program  in  well-defined  DATA  statements,  as  described 
in  Section  5.-2.  This  method  of  specifying  material  properties  iprovides 
flexibility  for  analyzing  various  candidate  materials  (at  the  expense  of 
program  recompilation  for  each  DATA  statement  alteration)  while 
reducing  the  number  of  inputs  required  for  normal1  design  analyses. 

The  case  wall  thickness  computed  from  Equation  3-1  is  compared 
with  a  specified  minimum  wall  thickness  based  on  fabrication  limita¬ 
tions,  and  the  larger  value  is  employed  in  further  calculations. 
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Following  an  initial  estimation  of  the  required  cylindrical  case 

liner  thickness,  t.  ,  the  outer  diameter  of  the  propellant  charge  is 
^  c 

computed  from 


D{  =  Dm  -  2  (Tm  +  T1(.)  .  (3-2) 

The  final  design  value  of  liner  thickness  is  determined  through  iteration 
using  an  empirical  relationship  which  is  described  later. 

Equation  3-2  begins  an  iteration  wherein  the  nozzle  expansion 
ratio,  cylindrical  case  liner  thickness,  and  propellant  geometrical 
characteristics  are  established.  An  attempt  is  first  made  to  size  the 
nozzle  for  optimum  expansion  such  that  the  nozzle  exit  pressure  is 
equivalent  to  the  input  ambient  pressure.  If  this  is  not  feasible  within 
the  constraints  of  the  input  maximum  nozzle  exit  diameter,  DGm>  and 
thrust-dictated  throat  area  requirement,  the  nozzle  design  which  most 
nearly  approaches  optimum  expansion  within  these  constraints  is 
established.  In  the  latter  situation,  the  nozzle  exit  pressure  is  always 
greater  than  ambient;  i.  e.  ,  the  nozzle  is  underexpanded.  As  a  first 
trial  estimate,  the  nozzle  exit  pressure,  Pe,  is  set  equal  to  the  :nput 
ambient  pressure  Pamb*  The  corresponding  nozzle  expansion .rs’io,  €  , 
is  then  determined  from  the  relationship 


€  = 


M4' fe>’ N' -  (*)  TT  -  '*• 


V  -1 


3) 


3-2 


The  thrust  coefficient.  Cp  is  computed  from 


Gi  = 


V  + 1 

Y-l 

f  2y  2 

(— 

/EeV  V  If 

Vv  +>, 

1  ! 

W  J  l 

(3-4) 

.  (Pe  ~  Pamb')'< 
+  - *c - 


The  divergence  loss  factor,  Xn,  is  defined  as 

Xn  =  j  (1  +  cos  a)  .  (3-5) 

For  analysis  of  contoured  nozzles  using  this  program,  an 
effective  exit  half-angle  of  17.  5  degrees  should  be  input,  as  described 
later.  The  nozzle  discharge  correction  factor,  q  is  usually  greater 
than  unity,  as  discussed  in  Reference  6,  but  may  range  from  0.  98  to 
1.15.  (A  value  of  1.  02  is  typical.  )  The  velocity  correction  factor, 
r|v,  ranges  between  0.  85  and  0.98,  with  an  average  near  0.  92  (Ref.  6). 

The  required  nozzle  throat  area  is  established  by 

At  =  F/(PC  Cf)  .  (3-6) 


The  corresponding  nozzle  exit  area,  Ag|,  is  computed  from 


The  required  nozzle  wall  thickness  at  the  exit  plane  (Tne)  *s  sized 
from  the  hoop  stress  relation 


Pe  Pei 


2<r 


n 


(3-8) 


The  variable  Fsn  is  a  program  input,  and  <rn  is  specified  within  the 
program  by  a  DATA  statement.  Again,  the  computed  wall  thickness 
from  Equation  3-8  is  compared  with  a  specified  minimum  thickness 
based  on  fabrication  constraints,  and  the  larger  value  is  selected  for 
further  calculations.  The  outside  diameter  at  the  nozzle  exit  plane  is 
then  computed  from 


Dep  =  Dei  +  2(Tne  +  Hx)  (3-9) 

where  Tix,  the  insulation  thickness  at  the  nozzle  exit  plane,  is 
computed  using  procedures  described1  later  in  this  section  (Equation  37). 
If  the  computed  outside  exit  diameter  as  less  than  the  input  maximum 
value,  the  design  analysis  continues  with  the  computed  value.  More¬ 
over,:  if  this  test  is  satisfied  on  the  first  iteration,  with  Pg  =  Pamb» 
the  nozzle  is  designed? for  optimum  expansion.  If,  however,  the  com¬ 
puted  outside  exit  diameter  exceeds  the  input  maximum  value,  the  exit 
pressure  estimate  is  increased  slightly  and  the  analysis  returns  to 
Equation  3-3.  The  iteration  continues  until  a  satisfactory  exit  diameter 
and  corresponding  exit  pressure  and  expansion  ratio  are  obtained. 


Following  the  establishment  of  the  nozzle  throat  and  exit  areas, 
the  interna!  ballistics  analysis  and  propellant  geometrical  characteriza¬ 
tion  proceed.  The  mass  discharge  coefficient  is  computed  from 


(3-10) 


The  independent  variables  y»  R>  Tc  are  program  inputs.  The  average 
mass  flow  rate,  m,  discharged  by  the  nozzle  is  computed  from 

-• 

m  =  CDPcAt  .  (3-11) 

This  discharged  mass  flowrate  must  equal  the  rate  of  mass  generation 
at  the  propellant  burning  surface,  which  is  determined  from 

m  =  rbAfcPp  (3-12) 

where  the  average  propellant  burning  rate  is  governed,  for  nonerosive 
flow  conditions,  by  the  well-known  relation 

rb  =  aPcn  .  (3-13) 


Equating  the  mass  generation  rate  (Equation  3-1 1 )  to  the  mass  discharge 
rate  (Equation  3-10)  yields,  after  rearrangement,  the  following  relation¬ 
ship  for  computing  the  required  average  burning  surface  area: 


Ab  = 


CDPcAt 

rbPp 


(3-14) 


The  initial  motor  port  area  is  established  by  the  input  port-to- 
throat  area  ratio  and  the  previously  computed  throat  area 

Ap  =  Apt  At  .  (3-15) 

The  initial  port-to-throat  area  ratio,  Apt,  must  be  input  as  unity  or 
greater  to  assure  a  nonerosive  flow  condition  inside  the  motor  port. 

For  tapered  grains,  Apt  is  the  length-averaged  value.  The  cross 
sectional  area,  Aw,  of  the  propellant  web  at  ignition  is  computed  from 
the  geometric  relationship 

Aw  =  (Af  -  Ap)  Fp  (3-16) 

where  Fp  is  the  fraction  of  propellant  burned  prior  to  web  burnthrough. 
Values  of  Fp  usually  range  from  0.  9  to  1. 0  (Refs.  2  and  3).  The 
corresponding  area,  As,  of  residual1  slivers  (active  or  inert)  remaining 
after  web  burnthrough  is  established  by  the  expression 


Ag  =  Af  -  A 


w 


(3-17) 


The  propellant  web  thickness  is  determined  from 


Tw  = 


Df  F. 


w 


(3-18) 


As  described  previously,  Fw  represents  the  input  propellant  web 
fraction.  Values  of  Fw  typically  range  from  0.  2  to  0.  5  for  star  or 
wagonwheel  designs  and  from  0.  28  to  0.  9  for  slottedrtube  designs. 

Star  and  wagonwheel  configurations  are  usually  employed  with  slow  or 
medium  burning  rate  propellants  where  relatively  small  burn  times  and 


large  thrust  levels  are  required.  Slotted  tube  configurations  yield 
large  motor  mass  fractions;  they  are  usually  employed  with  relatively 
slow  burning  propellants  for  long  burn  time,  low  thrust  applications  or 
with  fast  burning  propellants  for  short  burn  time,  high  thrust  applica¬ 
tions. 

The  motor  burn  time  is  computed  from 

tb  =  ^  •  (3-19) 

rb 

At  this  point  in  the  analysis,  the  required  liner  thickness  adjacent  to  the 
cylindrical  case  wall,  T|  ,  may  be  recalculated  more  accurately 
using  the  proven  empirical  expression  developed  by  Thiokol  Chemical 
Corporation  (Ref,  7). 

Tic  =  0.02  Dm0'5  (tb  +  l)0*1  (1  ^  Fp)0-25  •  (3-20) 

The  liner  thickness  computed  from  Equation  3-20  is  compared  with 
the  previously  estimated  value  used  in  Equation  3-2.  If  the  computed' 
liner  thickness  differs  by  more  than  a  specified  amount  from  the 
previous  estimate,  the  analysis  returns  to  Equation  3-2,  with  the 
computed  value  used  as  the  new  estimate.  This  iteration  continues 
until  the  computed  and  estimated: values  converge. 

The  next  step  in  the  analysis  is  the  determination  of  the  required 
propellant  volume  and  corresponding  motor  length.  Before  motor 
length  may  be  computed,  the  amount  of  propellant  must  be  evaluated 
which  may  be  loaded  into  the  available  forward' and  aft'head-end 
volumes,  if  head-end  webs  are  desirable.  (The  computer  program  is 
formulated:  such  that  the  analysis  of  head-end  webs  is  optional.  )  For 
cases  in  which  a  forward  head-end  web  analysis  is  requested,  a 


fractional  ellipsoidal  propellant  volume,  conforming  to  the  head 
closure  shape,  is  considered.  A  cylindrical  perforation  is  allowed 
for  protrusion  of  the  ignitor  through  the  head -end  web.  The  analysis 
of  the  propellant  in  the  propellant  volume  in  the  aft  head  is  similar  to 
that  for  the  forward  head,  except  the  cross-sectionalarea  of  the  center 
perforation  is  equivalent  to  the  port  area. 

Following  evaluation  of  the  combined  forward  and  aft  head-end 
propellant  volume,  vph'  and  corresponding  average  burning  surface, 
Ab^,  if  applicable,  the  required  average  burning  surface,  Abc,  within 
the  cylindrical  motor  length  is  computed  from 

Abc  =  Ab  -  Abh  -  (3-21) 

The  required  total  volume  of  the  propellant  charge,  Yg  ,  including 
residual  slivers  (active  or  inert),  in  the  cylindrical  portion  may  now 
be  computed  from  the  geometrical  relationship 


and  the  volume  of  residual  slivers,  Vs,  is  computed  from  the  geometri¬ 
cal  relationship 

VS  =  Vpsc  (1  -  Fp)  .  (3-25) 

Because  of  the  simple  propellant  geometry  selected  for  the  head-end 
web,  residual  slivers  there  will  be  negligible. 

The  required  cylindrical  length  of  the  motor,  Lc,  may  now  be 
computed  from 

Lc  =  VPcw/Aw  .  (3-26) 

The  propellant  and  residual  sliver  weights  are  evaluated  by  multiplying 
the  previously  computed  volumes  by  their  respective  densities.  For 
active  slivers,  the  propellant  density  is  used,  whereas  for  inert  slivers 
a  density  of  0.  02  lbm/:in3,  which  is  typical  of  conventional  phenolic/ 
microballoon  filler  material  is  used. 

3.  2  DESIGN  OF  INERT  COMPONENTS 

Sufficient  information  is  now  available  to  permit  computation  of 
weights  and  dimensions  of  the  remaining  constituent  inert  components 
of  the  rocket  motor.  These  calculations  are  performed  within- the 
computer  program  by  Subroutines  INERT,  NOZZ,  and  INS UL.  In  the 
interest  of  brevity,  only  the  final  formulations  of  the  pertinent  working 
equations  will  be  presented  herein.  The  reader  is  referred  to  the  cited 
references  for  the  appropriate  derivations.  The  nomenclature  for  the 
type  of  motor  case  considered  herein  is  defined  in  Figure  3—1.  Also, 
the  nomenclature  which  applies  to  the  conical  and  contoured  nozzle 
analyses  is  defined  in  Figures  3-2  and  3-3,  respectively. 


FIGURE  3-1.  MOTOR  CASE  NOMENCLATURE 


Nozzle  Structure 


FIGURE  3-2.  CONICAL  NOZZLE  NOMENCLATURE 


zle  Structure 


FIGURE  3-3.  CONTOURED  NOZZLE  NOMENCLATURE 


3.2.  1  Forward  Head  Membrane 

The  required;  thichness,  t^,  of  the  forward  head  membrane  is 
sized  by  the  resultant  of  the  meridional  and  tangential  stresses  as 
discussed  in  Reference  8.  Thus, 
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(3-27) 
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The  design  chamber  pressure,  Pp,  in  taken  as  1.5  times  the  average 
chamber  pressure.  The  value  of  Lh0  is  dictated  by  the  input  motor 
diameter  and  the  head-end  ellipse  ratio,  (3,  whichutypically  ranges 
from  1.  0  to  2.  0.  The  values  of  the  inside  major  axis,  D^,  and  semi- 
minor  axis,  L^,  are  determined  from  the  input  motor  diameter, 
ellipse  ratio,  and  membrane  thickness.  The  weight  of  the  forward  head 
membrane  structure  is  then  computed  from 

Wfhs  =  0*  167  Pmir  M%0  Dm2  -  D^-2 )  .  (3-28) 

3.  2.  2  Forward  Head  Insulation 

The  weight  of  the  forward  head  insulation  is  computed  from  the 
following  empirical  relationship  (Ref.  9): 


Whi  =  1.93  (TO-10)  Dm2 
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where 


* 


(3 -29a) 


3-13 


I 

3.  2.  3  Ignitor  Boss 

The  weight  of  the  ignitor  boss,  W.,  ,  on  the  forward  head  i.- 

lb 

estimated  from  the  empirical  expression  (Ref,  9): 
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(3-30) 


3.2.4  Ignitor 

The  weight  of  the  ignitor  ,  Wj[  ,  is  estimated  from  the  follow  - 
ing  empirical  relationship  which  was  derived  through  correlation  of 
reported  ignitor  weights  (Ref.  4)  for  previously  designed  and  fabricated 
rocket  motors: 


W. 

ign 
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.(3-31) 


Equation  3-31  was  found  to  produce  an  excellent  fit  of  reported  ignitor 
weight  data  for  rocket  motors  with  an  average  burning  surface  area 
ranging  between  3500  and  44000  in2. 

3.2.5  Cylindr i  cal  Ca s e  Str uc tur e 

The  weight  of  the  cylindrical  portion  of  the  motor  case  struc¬ 
ture,  Wcs>  is  computed  from  the'  relationship 


(3-32) 
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3.2.6  Cylindrical  Case  Liner 

The  weight  of  the  liner  adjacent  to  the  cylindrical  motor  c-.sc 
Wc/»  is  determined  from 
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(3-^3) 


For  the  computer  program  described  herein,  a  typical  liner  density, 
p  ,  of  0.  06  lbm/in3  is  employed 


3.2.7  Aft  Head  Membrane 

The  weight  of  the  aft  head  membrane  structure,  W^,  is 
estimated  from  the  following  empirical  relationship  (Ref.  9): 
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3.2.8  Aft  Head  Insulation 

The  weight  of  the  aft  head  insulation,  'W  is  scaled  as  1.54 
times  that  of  the  previously  computed  forward  head  insulation  as 
recommended  by  Reference  9. 


3.  2.9  Aft  Head  Nozzle  Attachment 

The  weight  of  the  aft  head  nozzle  attachment  boss,  Wab,  is 
estimated  using  the  following  relationship  (Ref.  9): 


W  .  =  10.26  P_  D  2  S 
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3.2.10  Nozzle 


Either  conical  or  contoured'nozzle  designs  may  be  analyzed 
using  this  program.  Conical  nozzles  of  the  type  illustrated  in  Figure 
3-2  are  designed  using  nondocumented  procedures  developed  by 
Mr.  A.  R.  Maykut  of  the  U.  S.  Army  Missile  Command  (MICOM)  at 
Redstone  Arsenal.  The  validity  of  these  procedures  has  been  proven  at 
Redstone  Arsenal  through  numerous  design  analyses  and  reproductions 
of  previous  operational  nozzle  designs.  Contoured  nozzle  length  and 
weight,  when  applicable,  are  scaled  from  a  corresponding  design  of 
a  conical  nozzle  with  equivalent  throat  area  and  fixed  exit  cone  half¬ 
angle  of  17;.  5  degrees  using  procedures  described  in  Reference  10.  A 
more  detailed  description  of  the  contoured  nozzle  scaling  procedure  is 
presented  later. 

Discussing  first  the  conical  nozzle  analysis,  the  entrance  wall 

thickness,  Te>  of  the  nozzle  at  the  aft  head  attachment  radius  is  compu¬ 
ted  from 


T  = 
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r  P_ 

op  D 

or  cos  <b~- 
n  T 


(3-36) 


The  parameter  r  usually  is  governed  by  the  size  of  the  aft  head 
opening  which  is  required  for  insertion  and  extraction  of  the  propellant 
mandrel.  For  this  computer  program,  the  nozzle  entrance  cone  half¬ 
angle,  <}> ,  is  fixed  at  60  degrees. 


The  thickness  computed  by  Equation  3-36  is  compared  with  a  specified 
minimum  value  which  is  dictated  by  fabrication  constraints,  and  the 
larger  of  the  two  values  is  used  in  further  calculations. 

The  next  step  in  the  nozzle  analysis;  is  the  determination  of  the 
required  insulation  thickness  at  the  nozzle  entrance.  The  procedure 
which  is  employed  in  sizing  the  insulation  thickness  is  described  in 
Reference  5.  A  thermally  thin  structural  wall  is  assumed  to  be  pro¬ 
tected  from  the  combustion  gases  by  a  refractory  insulation.  The 
temperature  of  the  exposed;  insulation  surface  is  assumed  to  be  equivalent 
to  that  of  the  adjacent  combustion  gases,  T^,  and  the  temperature  rise 
of  the  protected  structure  during  motor  operation  is  assumed  to  be 
much  smaller  than  T  .  The  predicted  temperature  rise  of  the  struc¬ 
tural  wall  during  time  t  is  computed  from  the  general  expression 
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where  0  is  the  time  constant  (t^2 /tt2  q^).  The  function 


in  Reference  5  as 


is  reported 
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Equation  3-37  cannot  be  solved  directly  for  because  of  the  implicit 
relationship  of  with  the  series  function  f(t/0).  Therefore,  an  itera¬ 
tive  solution  is  accomplished  in  which  Equation  3-37  is  solved  repeti¬ 
tively  to  compute  values  of  predicted  temperature  rise  for  progressively 
improved  estimates  of  r^.  The  least  value  of  for  which  the  predicted 
structural  temperature  rise  is  maintained  below  a  specified1  limit  of 
10°R  is  selected  for  design  application. 
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The  first  20  terms  are  utilized  in  the  series  evaluation  of  f(t/0). 


Ablation  of  the  insulation  is  not  accounted  for  in  the  determination  of 
the  required  insulation  thickness.  Therefore,  the  selected  thickness 
should  be  conservative.  An  approximate,  but  nonconservative,  means 
of  accounting  for  ablation  would  be  to  substitute  the  insulation  ablation 
temperature  for  T^  in  Equation  3-37:. 

The  weight  of  insulation  in  the  nozzle  entrance  cone,  Wej.,  is  com¬ 
puted  from 
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The  weight  of  the  nozzle  entrance  cone  structural  material,  Wes,  is 


computed  from 
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The  weight  of  the  nozzle  throat  insulation  is  computed  from  the  relation 
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(cos  4>  +  cos  a  )  +  r^t  ^k^r^  .  (3-41) 


The  thickness  of  the  throat  structural  shelf,  Tts,  is  sized  from  the 
following  expression  which  accounts  for  the  reductioninpressure  load¬ 
ing  of  the  structural  shell  due  to  the  load  carried  by  the  throat  insert 


material: 
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The  design  pressure  Pj}  in  Equation  3-42  should:  rigorously  be  based 
upon  a  gas  pressure  whose  value  lies  between  the  expected  values  at 
the  throat  and;  in  the  combustion  chamber,  since  the  inner  surface  of 
the  throat  insert  is  exposed  to  pressures  in  this  range.  However,  for 
design  conservatism,  the  value  of  Pj)  used  in  this  program  is  computed 
as  1.5  times  the  chamber  pressure.  The  corresponding  weight  of  the 
throat  structural  shell,  Wts«  is  determined  from 
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From  geometrical  considerations,  the  weight  of  the  throat 
insert,  Wj,  is  established- as 
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The  thickness  of  the  exit  cone  structure  at  the  downstream  end  of  the 


nozzle  throat  insert, rxs,  is  sized  from 
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(3-45) 


The  design  pressure  Pj}x  4s  defined  as  1.5  times  the  local  pressure  in 
the  nozzle  exit  cone.  The  structural  thickness  computed  from  Equa¬ 
tion  3-45  is  compared  with  a  specified  minimum,  value  which  is  dictated 
by  fabrication  constraints,,  and  the  larger  of  the  two  values  is  selected 
for  design  application.  The  corresponding  weight  of  the  exit  cone 
structure,  Wxs,  is  computed  from 
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The  weight  of  the  nozzle  exit  cone  insulation,  W  is  determined  from 
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The  weight  of  the  attachment  boss,  VI  ,  which  is  an  integral  part  of 

nb 

the  nozzle  structure,  is  estimated  from 
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From  the  previously  computed  nozzle  component  weights,  the  total 
nozzle  weight,  W^,  is  computed  as  the  sum 


W  =  W  .  +  W  +  W\.  +  W  +  W  +  W  +  W  .  4  W  ,  .  (3-49) 
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The  nozzle  length,  L  ,  for  conical  entrance  and  exit  geometries, 


n 


may  be  computed  from 
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Throat 


Exit 


Vf'op-ro|'c,n't’  +  rtkt<sM>+  4ina)x  -  r^jctr.o  .  (3-50) 


Equations  3-36  through  3-50  were  derived  for  nozzles  with 
conical  entrance  and  exit  geometries. 
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However,  most  modern-day  rocket  motors  which  are  designed  ,for 
flight  operation  employ  nonconical  contours  which  are  optimized  for 
improved  performance  and  reduced  length.  Because  of  the  difficulty 
associated  with  detailed  design  of  a  contoured  nozzle,  the  simplified 
empirical  procedure  described  in  Reference  10  for  approximating 
contoured  nozzle  weight  and  length  from  a  corresponding  conical 
nozzle  design  was  employed  in  this  program. 

To  calculate  contoured  nozzle  weight  using  t’>e  previously 
described  approximate  method,  the  following  assumptions  are  made: 

•  The  weights  of  conical  and  contoured  nozzles  are 
identical  from  the  forward  attachment  boss  to  an 
expansion;  ratio  downstream  of  the  throat  of  3.  0  for 
nozzles  of  equal  throat  area.  (Nozzle  shape  within 
this  region  is  relatively  independent  of  downstream 
shape. ) 

•  The  nozzle  weights  downstream  of  an  exit  expansion 
ratio  of  3.  0  are  identical  for  contoured  and  conical 
nozzles  of  equal  surface  area, 

The  two  types  of  nozzles  are  assumed  to  have  identical  throat  areas, 
chamber  presures,  burning  durations,  and  propellant  flame  tempera¬ 
ture.  The  family  of  contoured  nozzles  for  which  this  procedure  applies 
is  described  by  O.  J.  Demuth  in  Reference  11, 

To  determine  the  conical  nozzle  surface  area  which  is  equiva¬ 
lent  to  that  of  a  contoured  nozzle  having  the  same  throat  area,  an 
equivalent  conical  nozzle  expansion  ratio,  «e,  is  derived  from  previ¬ 
ously  computed;  contoured  nozzle  data.  This  equivalent  expansion 
ratio  is  defined  as  that  of  a  conical  nozzle  with  a  17.  5-degree  exit 
cone  half-angle  which  has  the  same  surface  area  as  a  particular  con¬ 
toured  nozzle  of  the  same  throat  area.  Contoured  nozzle  weights  are 
then  computed  from  the  previously  described  conical  nozzle  equations 
(Equations  3-36  through  3-50)  when  the  applied  value  of  the  inside  radius 


of  the  throat  exit  plane,  re,  is  based  on  the  equivalent  expansion  ratio, 
Cg,  and  the  input  nozzle  exit  cone  half-angle  is  taken  as  17.  5  degrees. 
The  curve  of  €e  as  a  function  of  «  (Ref.  10),  which  is  incorporated 
into  this  computer  program,  is  presented  in  Figure  3-4.  Figure  3-4 
applies  for  a  Demuth-type  contoured  nozzle  with  an  initial  divergence 
half-angle  (immediately  downstream  of  the  throat)  of  32.5  degrees 
and  an  expansion  ratio  at  the  point  of  parallel  flow  of  30.  A  review  of 
previously  designed  contoured  nozzles  (Ref.  4)  revealed  that  this 
particular  exit  shape  is  generally  applicable  for  relatively  large  tactic- 
cal  rocket  motors.  Of  course,  flight  nozzles  are  truncated  at  expan¬ 
sion  ratios  slightly  less  than  that  required  for  parallel  flow  with  little 
performance  penalty.  Additionally,  Figure  3-4  is  based  on  a  ratio  of 
combustion  gas  specific  heats  of  1.  2  and  a  ratio  of  throat  radius  of 
curvature  to  throat  radius  of  1.  2. 

Direct  calculation  of  the  length  of  a  contoured  nozzle  is  again 
quite  difficult.  To  simplify  the  calculations  in  this  program,  contoured 
nozzle  length  is  related  empirically  to  a  17.  5-degree  conical  nozzle  of 
the  same  throat  area  and  expansion  ratio.  In  determining  contoured 
nozzle  length,  the  length  from  the  throat  to  the  exit  plane  is  first  cal¬ 
culated  as  for  a  17.  5-degree  conical  nozzle.  This  length  is  then  multi¬ 
plied  by  an  empirically  derived  factor,  Lcont/Lcon,  for  the  appropriate 
contoured  nozzle  expansion  ratio.  The  curve  of  Lconj./Lcon  as  a  func¬ 
tion  of  €  (Ref.  10),  which  is  employed  in  the  computer  program  de¬ 
scribed  herein,  is  presented  in  Figure  3-5,  The  curve  shown  in 
Figure  3-5  is  for  the  same  Demuth-type  nozzle  shape  as  discussed 
previously  for  nozzle  weight  calculations.  The  corrected  length  between 
the  throat  and  the  nozzle  exit  plane  is  then  added  to  the  entrance  length, 
between  the  forward  attachment  boss  and  the  throat,  to  determine  the 
total  contoured  nozzle  length. 


Equivalent  Conical  Nozzle  Expansion  Ratio 


FIGURE  3-4.  EQUIVALENT  CONICAL  EXPANSION  RATIO  AS  A  FUNCTION 
OF  CONTOURED  NOZZLE  EXPANSION  RATIO 
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Contoured-Nozzle  Expansion  Ratio,  e 


Lcont/Lcon  _AS  A  FUNCTI0N  0F  CONTOURED  NOZZLE 
EXPANSION  RATIO 


* 


The  previously  described  procedure  could  be  extended  for 

analyzing  other  Demuth-type  exit  shapes,  with  different  values  of 

«.  and  «  for  parallel  flow,  by  incorporating  additional  curves  of  ce 

and  L  /  L  as  a  function  of  (  from  Reference  10. 
cont/  con 


3.2.  11  Motor  Attachment  Skirts 

The  weight  of  the  forward  attachment  skirt,  WfSl  is  estimated 
from  the  following  relationship  recommended  by  Reference  9: 


W.  =  C,  D  2 
fs  1  m 


g  w 

maxi  PL 


).  215  (L  +D  /P\ 

\  c  m  )  . 


(3-51) 


The  maximum  longitudinal  acceleration  of  the  stage  under  considera¬ 


tion,  gm.v. ,  is  assumed  to  occur  at  stage  burnout: 


g  «  F/W.  . 
max.  bo 

l 


(3-52) 


Now,  at  this  stage  of  the  analysis,  is  unknown,  since  its  value 

depends  upon  the  undetermined  weights  of  the  attachment  skirts  and 

interstage  structure.  However,  these  weights  are  usually  relatively 

small  in  comparison  with  the  remaining  inert  component  weights. 

Therefore,  an  iterative  solution  is  accomplished  in  which  g  is 

"max. 

x 

first  estimated  using  a  value  of  which  includes  only  the  previously 
computed  weights  of  the  motor  case,  nozzle,  ignitor,  liner,  insula¬ 
tion,  and  the  stage  payload.  This  value  of  g  is  utilized  in 

rna.x^ 

the  calculation  of  the  attachment  skirt  and  interstage  weights.  The 

value  of  g  is  then  recalculated  using  the;estimated  skirt  and  inter- 
°max. 
n 

stage  weights  included  in  This  iteration  is  continued  until  the 


3-25 


computed  values  of  Wfco  '(°r  gmaxi)  on  successive  iterations  are 
equivalent  within  a  specified  tolerance.  At  this  point,  the  correct  values 
of  the  attachment  skirt  and  interstage  weights  are  considered  to  be 
established. 

The  weight  of  the  aft  attachment  skirt,  Was,  is  estimated  from 


was  =  o.  0055  Dm! 


Smax;  i  W0 


'0.215  (lc  +  l)0-5 

—  DZ - —  +  1  * 


(3-53) 


In  the  evaluation  of  Was,  the  value  of  gmax^  j  must  be  estimated  from 
a  separate  analysis  of  stage  i-1,  and  WQis  estimated  as  part  of  the 
previously  described  iteration  oh  Wj^q. 

3.2.12  Interstage  Structure 

The  interstage  structural  weight,  W int»  Is  estimated,  when  re¬ 
quired,  from  the  following  relationship  which  is  recommended  by 
■Reference  9: 


Wint  =  0-  15?  Dm  +  L„[+1  +  Lint) 


gmaxi  WPL  f 0.  21 5 


h :  M  +  if 

Dm  J) 


(3-54) 


The  evaluation  of  gmaxj  and  Wint  is  included  in  the  previously  described 
iteration  on  W_b0. 


3.  3  ROCKET  MOTOR  PERFORMANCE 
CHARACTERIZATION 

Sufficient  preliminary  design  information  has  now  been  computed 
from  the  previously  described  analysis  to  permit  prediction  of  rocket 
motor  performance.  After  setting  the  sum  of  the  previously  computed 
inert  component  weights  equal  to  Wjp,  the  gross  stage  weight,  W5,  is 
computed  from 


(3-58) 


Stage  mass  fraction,  |x,  is  defined  as 


Wp 

11  ■  Wo  -  W:PL  ‘ 

Rocket  motor  mass  fraction,  X,  is  defined  as 


X  = 


Ie _ 

Wint  “  WPL  * 


(3-59) 


(3-60) 


The  delivered  specific  impulse,  ISp^,  computed  from 

IsPd =  c^;  •  (3_61) 

The  predicted  ideal  velocity  increment  produced  by  the  stage,  neglect¬ 
ing  drag  and  gravity  losses,  is  computed  from 

AV:  -  ISpd  gc  in  Ea.  '  <3-62) 

Thrust-to-weight  ratios  at  ignition,  Aign,  and  burnout,  Abo,  are  com¬ 
puted  from  Equations  3-63  and  3-64,  respectively: 


Aign  =  F/Wq 

(3-63) 

Ajjo  =  F/Wbo  • 

(3-64) 

This  completes  the  description  of  the  rocket  motor  design 
analysis  as  it  is  programmed  for  the  digital  computer.  The  program 
is  formulated  in  such  a  manner  that  parametric  tradeoff  .analyse  s  of 
pertinent  design  variables  may  be  readily  accomplished  by  stacking 
input  cases  back-to-back.  Moreover,  an  automatic  variable  adjusting 
procedure,  which  would  determine  the  required  value  of  a  specific 


independent  variables  (e.  g. ,  motor  diameter,  chamber  pressure, 
propellant  web  fraction,  etc.  )  to  produce  a  desired  motor  performance 
or  design  characteristic  (e.g.,  mass  fraction,  motor  length,  or 
velocity  increment)  may  be  easily  incorporated  into  the  program  if 


4.  PROGRAM  APPLICATION 


4.  1  METHODOLOGY 

The  computer  program  described  herein  may  be  utilized  to 
generate  both  point  and  optimized  preliminary  designs  of  solid  propellant 
rocket  motors  through  the  proper  selection  and  manipulation  of  appropri¬ 
ate  input  variables.  Typical  design  problems  and  corresponding  methods 
for  using  the  program  to  solve  these  problems  are  presented  in  Table  4-1. 
The  computer  program  contains  no  automated  variable  optimization 
scheme.  Therefore,  the  optimization  of  pertinent  design  variables  must 
be  accomplished  using  "brute  force"  techniques  wherein  graphs  are 
plotted  of  computed  performance  and/or  weight  characteristics  which 
result  from  manual  perturbation  of  the  variables  of  interest.  Optimum 
values  of  design  variables  which  yields  either  maximum  performance  or 
minimum  weight  characteristics  are  then  selected  from  these  graphs. 

For  example,  Problem  1  of  Table  4-1  illustrates  the  technique  which 
would:  be  employed  to  establish  the  optimum  chamber  pressure  for  a 
minimum  weight  design,  assuming  all  other  input  variables  are  fixed. 

Problem  2a  of  Table  4-1  illustrates  how  the  designer  would 
determine  the  optimum  motor  diameter  and  corresponding  chamber 
pressure  which  would  be  required  to  provide  a  given  motor  length. 
Similarly,  Problem  2b  illustrates  the  technique  which  would  be  used  to 
determine  the  optimum  motor  length  and  corresponding  chamber  pressure 
which  would  yield  a  given  motor  diameter.  The  latter  two  examples  are 
frequently  encountered  in  the  design  of  tactical  and  strategic  missiles 
which  must  be  contained  within  fixed  envelope  dimensions. 

Problem  3  illustrates  how  the  program  would  be  employed  to 
evaluate  the  variation  of  stage  impulsive  velocity  with  motor  diameter, 
expansion  ratio  and  payload  weight. 


TABLE  4-1.  PROGRAM  APPLICATION  TO  TYPICAL  DESIGN  PROBLEMS 
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Problem  4  illustrates  how  the  sizingtprogram  described  herein 
may  be  used  in  conjunction  with  a  stage  optimization  routine  to  develop 
a  near-optimum  preliminary  design  of  the  propulsion  components  of  a 
multistage  vehicle.  A  typical  stage  optimization  routine,  which  was 
computerized  by  Mr.  J;.  H.  Dobkins  of  Teledyne  Brown  Engineering,  is 
described  in  Appendix  B.  Results  from  this  routine,  which  is  entitled 
the  Missile  Optimization  Program  (MOP),  have  been  demonstrated  to 
agree  closely  with  those  from  other  more  elaborate  analyses.  Since 
MOP  requires  ideal  impulsive  velocity  as  input,  the  designer  must 
estimate  velocity  losses  which  will  result  from  drag  and  gravitational 
effects.  Additionally,  input  values  of  mass  fraction  and  specific  impulse 
for  each  stage  must  be  specified.  Minimized  gross  vehicle  weight  and 
optimized  propellant  weight  distributions  for  a  prescribed  number  of 
stages  are  then  computed  using  MOP  for  the  ranges  of  payload  weights 
and  ideal  impulsive  velocities  of  interest.  For  selected  point  values  of 
ideal  impulsive  velocity  and  payload  weight,  the  designer  computes 
corresponding  values  of  required' total  impulse  and  thrust  (for  a  pre¬ 
scribed  burn  time)  for  each  stage.  The  sizing  program  may  then  be 
utilized  to  establish  optimum  diameters  (or  lengths),  chamber  pressures 
and  expansion  ratios  for  all  of  the  stages  as  discussed! previously. 

Following  the  preliminary  design  of  a  single  or  multistage  /ehicle, 
the  evaluation  of  its  predicted  flight  performance  is  usually  desired. 

The  trajectory  analysis  which  is  -required  to  evaluate  flight  performance 
is  beyond  the  scope  of  this  investigation.  However,  vehicle  design 
characteristics,  e.  g.r  weights,  mass  flowrates,  thrust  levels, 
envelope  dimensions,  which  are  desired  from  the  sizing  program  may 
be  input  into  a  trajectory  simulation  computer  program  to  predict 
flight  performance.  The  trajectory  simulation  program  described  in 
Reference  12  is  typical  of  many  applicable  programs  which  are  available. 


4.  2 


SAMPLE  PROBLEM 


To  demonstrate  the  validity  of  the  computer  program,  a  sample 
problem  was  processed  which  consisted  of  an  analytical  reproduction  of 
the  Thiokol  TX-354-3  (Castor  II-A)  rocket  motor  design  described  in 
Reference  13.  This  motor  was  chosen  for  discussion  herein  because 
of  the  ready  availablity  of  unclassified  documentation  describing  its 
desigmand  performance  characteristics  in  detail.  Other  rocket  motor 
designs  described  in  Reference  4,  e.  g.  ,  Spartan,  were  also  analyzed 
successfully  using  this  program,  but  the  classified  nature  of  their  design 
characteristics  precluded  their  discussion  herein. 

The  cross-sectional  arrangement  of  the  propellant  grain  in  the 
TX-354-3  rocket  motor  is  illustrated  in  Figure  4-1.  The  grain  design 
is  basically  a  cylindrical  port  type  which  is  segmented  by  two  radial 
slots,  as  shown  in  Figure  4-1.  A  propellant  web  fraction,  Fw,  of  0.  76  was 
deduced  from  the  published  geometric  characteristics  of  the  propellant 
charge.  The  TX-354-3  nozzle  design  (Figure  4-2)  utilizes  a  conical 
expansion  geometry,  a  graphite  throat  insert,  and  requires  no  entrance 
structure.  The  nozzle  is  designed  for  operation  at  near-vaccum  condi¬ 
tions.  AISI  4130  steel  is  utilized  as  the  structural  material  for  the 
rocket  motor  chamber  and  nozzle  body. 

The  required  input  parameters  corresponding  to  the  TX-354-3 
rocket  motor  were  translated  into  a  set  of  motor  design  arid  performance 
characteristics  using  the  computer  program  described  herein.  The 
computer  printout  of  the  input  and;  output  parameters  of  the  sample 
problem  analysis  is  presented  in  Figure  4-3.  For  comparison,  published 
values  (Ref.  13)  of  pertinent  design  and  performance  parameters  of 
the  TX-354-3  motor  are  included1  in  parentheses  adjacent  to  the  corre¬ 
sponding  computer  values.  Comparison  of  the  published  and  computed: 
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FIGURE  4-3.  SAMPLE  PROBLEM  SOLUTION 


values  reveals  certain  discrepancies  which  must  be  discussed.  For 
example,  the  computer  nozzle  weight  of  170.5  lbm  is  much  less  than 
the  published  value  of  539  lbm.  This  large  discrepancy  is  difficult  to 
explain  without  a  knowledge  of  the  details  of  the  original  structural 
analysis  from  which  the  nozzle  was  designed.  However,  it  may  be 
speculated  that  bending  stresses  from  predicted  lateral  flight  loads 
could  ha.ve  dictated  larger 'required  structural  thicknesses  than  those 
which  resulted  from  the  pure  hoop  stress  equations  employed  in  this 
computer  program.  (The  credibility  of  the  nozzle  design  procedures 
employed  herein  was  demonstrated  through  the  almost  exact  analytical 
reproduction  of  the  existing  Spartan  booster  nozzle.;):  The  only  other 
conceivable  explanation  for  this  discrepancy  is  that  relatively  large 
design  safety  factors  (»  1.5)  could  have  been  employed  in  the  original 
design  analysis  of  the  TX-354-3  nozzle. 

The  only  other  significant  discrepancies  between  published  and 
computed  motor  design  parameters  are  in  total  motor  length,  total 
inert  parts  weight,  and  mass  fraction.  The  discrepancy  between  the 
published  motor  length  of  245  inches  and  the  computed  value  of  242.  6  inches 
results  primarily  from  the  neglect  in  the  computer  program  of  the  reported 
lengths  associated  with  the  pyrogen  and  nozzle  attachment  bosses  (Figure 
4-1).  The  discrepancy  between  the  published  total  inert  parts  weight  of 
1523  lbm  and  the  computed  value  of  811.4  lbm  results  primarily  from 
the  nozzle  weight  deficit  (=368.  5  lbm),  the  chamber  weight  deficit  (92  lbm), 
and  the  weight  deficit  caused  by  neglect  of  the  additional  insulation  weight 
required  by  the  two  radial  slots  in  the  propellant.  Additionally,  the 
discrepancy  between  the  published1  mass  fraction  of  0.  84  and  the  computed 
value  of  0.91  results  directly  from  the  previously  described  discrepancy 
in  total  inert  parts  weight. 


Additional  test  cases,  covering  a  variety  of  motor  configurations, 
should  be  processed  in  the  future  to  clearly  establish  the  accuracy  and 
general  applicability  of  the  computer  program. 
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5.  PROGRAM  OPERATING  INSTRUCTIONS 


The  previously  described  design  equations  were  coded;  into  a 
FORTRAN  computer  program  consisting  of  approximately  530  cards. 
Multiple  cases  may  be  submitted  simultaneously  by  stacking  the  input 
cards  back-to-back.  Only  the  input. parameters  which  change  need  be 
added  to  complete  each  new  case.  The  remaining  input  parameters 
always  revert  to  the  values  used  in  the  preceding  case. 


5.  1  INPUT  DATA 

All  program  input  variables  are  read  into  the  array  A(i),  i  = 

1,  27.  A  card-by-card  description  of  the  inputs,  including  symbol 
definition  and  format,  is  presented  below. 

•  Card  1,  Format  (12) 

A  Code  word  KTR 

A  KTR  is  the  number  of  input  variables  (one  per  card)  to 
be  read  following  this  card.  At  the  end  of  the  last  case, 
input  KTR  as  99  to  end  the  job. 

•  Cards  2  through  28,  Format  (12,  E15.8) 

A  These  27  cards  contain  the  27  input  variables  which  are 
read  into  the  array  A(i).  On  each  of  the  27  cards,  the 
proper  index,  i,  is  punched  into  columns  1  and  2  ,  and 
the  corresponding  element  A (i )  is  punched  into  columns 
3  through  17,  using  the  format  given  above. 

A  The  27  inputs,  A(i),  i  =  1,  27,  and  the  equivalent 
program  variables  are  defined  in  Table  5-1. 

•  Card  29,  Format  (12) 

A  Code  word  KTR 

A  This  is  the  same  code  word  as  described  above.  If 
another  case  is  to  follow,  KTR  is  the  number  of  input 
elements  A(i)  to  be  changed  from  the  values  used  in 
the  preceding  case,  if  no  more  cases  follow,  KTR 
should  be  punched  as  99- 


I 


TABLE  5-1.  INPUT  VARIABLE  DEFINITION 
FOR  CARDS  2  THROUGH* 28 


- 

ELEMENT 

EQUIVALENT 

PROGRAM 

VARIABLE 

DEFINITION 

UNITS 

* 

A(l) 

DM 

Rocket  motor  outside  diameter 

in. 

1 

A(2) 

PC 

'Chamber  pressure 

psia 

A(3) 

F 

Required  thrust 

F 

i 

1 

A(4) 

FW 

Pronellant  web  fraction 
(Ratio  of  web  thickness  to 
propellant  outside  radius) 

- 

A(  5) 

APT 

Port-to-throat  area  ratio 

— 

i 

A(6) 

PAM 

Ambient  pressure 

psia 

A(7) 

TC 

Propellant  flame  temperature 

°R 

r* 

A(8) 

RHOP 

Propellant  density 

lbm/in3 

L 

A(9) 

GAM 

Combustion  gas  ratio  of  specific 
heats 

-- 

L 

A(10) 

ATB 

Pressure  coefficient  in  pro¬ 
pel  lant.burning  rate  law 

— 

F 

A(ll) 

AMW 

Molecular  weight  of  combustion 
gas 

Ibm/mole 

A(12) 

AN 

Pressure  exponent  in  propellant 
burning  rate  law 

-- 

[ 

A(13) 

ETAD 

Mass  discharge  correction  factor 
for  nonideal  flow 

— 

I 

A(14) 

ETAV 

Exit  velocity  correction  factor 
for  nonideal  flow 

— 

* 

A(15) 

ALF 

Nozzle  exit  cone  half-angle  (for 

deg 

contoured  nozzle  analysis,  input 
effective  value  of  ALF  of  17.5 
degrees) 


l 


5-2 


TABLE  5-1.  INPUT  VARIABLE  DEFINITION 
FOR  CARDS  2  THROUGH  28  (Continued) 


ELEMENT 

EQUIVALENT 

PROGRAM 

VARIABLE 

DEFINITION: 

UNITS 

A(16) 

DEM 

Maximum  permissible  nozzle  exit 
"plane  diameter 

in. 

AO  7) 

FP 

Fraction  of  propellant  burned 
before  web  burn  through 

— 

’ 

A(18) 

SFM 

Motor  case  design  safety  factor 

— 

i 

A(19) 

SFN 

Nozzle  design  safety  factor 

— 

l 

A(20) 

WPL 

Stage 'payload  weight  (includes 
terminal  payload  and  upper 
stage  weight,  where  applicable) 

Ibm 

i 

A(21) 

WEX 

Miscellaneous  weights  to  be 
added  to  stage  weight  (e.g., 
thrust  vector  control  hardware) 

Ibm 

1 

l 

A(22) 

2LN 

Nozzle  length  of  next  higher 
stage 

in. 

j 

A(23) 

ZLI 

Interstage  clearance  between 
forward  head  of  stage  being 
designed  and  nozzle  exit  plane 
of  next  higher  stage 

in. 

f 

A(24) 

CODE 

Code  designating  which  stage 
i s  being  designed  (CODE  =  2.0 
for  terminal  stage  of  multiple 
vehicle  or  for  single-stage 
vehicle;  CODE  =1.0  for 
remaining  stages) 

* 

A(25) 

FLAG 

Set  FLAG  =0.0  when  head  end 
web(s)  are  not  present;  FLAG  = 

1.0  forward~Hiad  web  only; 

FLAG  =  2.0  for  forward  and  aft 

head  webs;  FLAG  =3.0  for  aft 
head  web  only. 


I 

I 
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TABLE  5-1.  INPUT  VARIABLE  DEFINITION 
FOR  CARDS  2  THROUGH  28  (Concluded) 


ELEMENT 

A(26) 

A(27) 


i 

} 


EQUIVALENT 

PROGRAM 

VARIABLE  DEFINITION 


NOZ  .Set  NOZ  =0.0  for  conical 

nozzle  analysis;  NOZ  =  1.0 
for  contoured  nozzle  analysis 

ISLIV  For  cases  where  FP  <1.0,  set 

ISLIV  =0.0  for  active  slivers; 
ISLIV  >0.0  for  inert  slivers 


5-4 


UNITS 


In  preparation  of  the  program  inputs,  FORTRAN  coding  sheets 
are  usually  completed  to  serve  as  a  guide  in  punching  the  appropriate 
values  onto  input  cards.  Example  input  coding  sheets,  including  input 
variable  definitions  and  typical  ranges  and/or  recommended  values  of 
the  inputs,  are  presented  in  Figure  5-1. 


FIGURE  5-1.  SAMPLE  INPUT  CODING  SHEET 


SAMPLE  INPUTS* 

PROGRAM 

VARIABLES** 

TYPICAL  RANGES  OR 

RECOMMENDED  VALUES 

f  s « 

7  10  15  20. 

2.7 

...  1  ....  1  ....  1 

_ 

ktr 

flumDer  of  inputs  to  follow  or 
gg  to  onrf  jnh 

/  +  .  3  / 

.  .  .  i  .  £+1 02.  .  .  1 

- 

DM 

Constrained  by  envelope 
limitations 

^2.+. '.(,2 

9  i  f  y-io  3  i 

PC 

300  to  2000  DSia 

-3  +  *  L  / 

7i>0  I  .  .  E.+10.S.  .  ,  i 

F 

Dictated  by  total  impulse, 
hum  time  requirements 

V.+--  76 

i  e  .-ho  o  .  i 

FU 

Sfetf.SWWi  SW8  °-5 

.5.+  *  1  2 

7.  .  .  I  .  .  .£.+10./.  .  .1 

APT 

1.1  to  3.5 

^.+••.5 

.  .  .  1  .  .  .£.+10.0.  .  .1 

PAM 

0.0  to  14.7  psia 

r.-hiAi 

V.2  l  .£.+10  */  .  I 

TC 

4000  to  6600  °R 

$±.•01 

.  1  .  .  .£.+10.0.  .  .  1 

RHOP 

0.053  to  0.07  lbm/in3 

9+  •  /  1 

£  i  £+10  /  .  i 

GAM  ._ 

1.12  to  1.2  _  . 

To*  \OS 

„  ,  £+10  0  , 

ATB 

0.002  to  0.1 

1  T+  •  23 

.  .  1  .  .  .£.+10.2.  .  ,  1 

. AMW 

25  to  29  . 

I.Z.t-27 

'  .  1  .  .  .£.+10.0.  .  .  1 

AN 

0.1  to  0.9 

1HSR 

iMHHHUHM 

ETAD 

0.98  to  1.15  (1.02  typ.) 

BPSR 

iimi nitiuu  jbwujj 

ETAV 

SHBB 

ALF 

WRMI 

IHHHHR99H 

DFM  — 

Constrained  by  envelope, 
expansion  ratio  nnntatiOiis 

Rffi 

FP 

0.9  to  1.0 

SFM 

1.5 

/  9+  *./|a 

i{  .  .  .  i  .  .  .£.+10./.  .  ,  i 

SFN 

1.5 

f 

20  +  -  a;; 

;/.0  i  EmoA.  .  .  i 

WPL 

Dictated  bv  mission  objectives 

! 

2  i  +  •  d 

.  .  i  £+100  1 

WEX  — .  - 

Additional  miscellaneous  weight 
(Fins.  TVC  Hardware,  etc.) 

-  -  i 

Z£.+  uCr 

_ !  ..  .  .£+10.0  .  1 

ZLN 

K  <}  iLk! 

gJ.tj.jg. 

fed 

i  .1.  .  E.+I0.0  .  .  L. 

zu  z 

3  -  10  in. 

1 

1  _ t  .£.+10. 1 .  .  .  i 

1 

★* 

J_. _ _  I  £+10./.  .  .1 

■3H 

** 

NOZ 

** 

in; 

ISLIV 

** 

m 

HHHHMMI 

KTR  _ 

End  of  job 

*  Values  listed  are  inputs  for  Sample  Problem  (Figure  4-3) 
** "Refer  to  page  5-1  and  Table  5-1  for  Definitions 


l* 
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5.  2  STORED  DATA 

Most  of  the  material  physical  properties  which  are  required  in 
the  motor  design  are  defined  through  the  use  of  FORTRAN  DATA 
statements  within  the  main  segment  of  the  program.  This  procedure 
was  adopted  to  reduce  the  number  of  inputs  which  would  be  required 
in  parametric  analyses  once  the  motor  case,  insulation,  and  nozzle 
materials  have  been  selected  for  a  particular  rocket  motor  design. 
Program  flexibility  is  not  impaired  significantly,  because  various 
structural  materials  may  be  readily  analyzed  by  simply  modifying 
the  appropriate  property  values  in  the  previously  mentioned  DATA 
statements  before  program  compilation.  These  property  values  are 
then  employed  by  the  program  until  further  changes  in  the  DATA 
statements  are  made. 

The  program  variables  which  represent  the  various  physical 
properties  are  well'defined  in  the  main  program  listing  in  the  Appendix 
by  several  COMMENT  cards.  Program  variable  definitions  and 
corresponding  typical  property  values  are  presented  in  Table  5-2, 

The  values  shown  in  Table  5-2  are  those  incorporated  in  the  program 
deck  which  is  listed  in  the  Appendix. 
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TABLE  5-2.  TYPICAL  MATERIAL  PHYSICAL  PROPERTIES 
AND  EQUIVALENT  PROGRAM  VARIABLES 


PROGRAM 

VARIABLE 

DEFINITION 

TYPICAL  VALUE 

r  RHOM 

Motor  case  density 

0.283  lbm/in3 

Motor1 

Case  { 

SIGM 

Allowable  tensile 
stress 

150,000  psi 

EC 

Young's  modulus  of 

30  (106)  psi 

case 

El 

Young's  modulus  of 
interstage  structure 

30  (106)  psi 

BETA 

Forward  and  aft  head 

1.0  -  1.6 

V. 

ellipse  ratio 

^RHOS 

Density 

0.283  lbm/in3 

SIGN 

Allowable  tensile 

150,000  psi 

Nozzle1  , 
Structure  ' 

stress 

;  ESUBS 

Young's  modulus 

30  (10e)  psi 

MUS 

Poisson's  ratio 

0.3 

^  CSUB2 

Specific  heat 

0.11  Btu/lbm-°R 

f  RHQI 

Density 

0.062  lbm/in3 

Nozzle2  j 
Insulation^ 

KSUB1 

Thermal  conductivity 

0.5  (10**5)  Btu/in. 

1  CSUB1 
^RHOT 

Specific  heat 

0.21  Btu/lbm-°R 

Density 

0.062  lbm/in3 

Throat3  ^ 
Insert 

ESUBT 

Young's  Modulus 

1.7  (106)  psi 

MUT 

Poisson's  ratio 

0.12 

SIGT 

Allowable  tensile 

3,355  psi 

stress 

Insert4  / 
Slivers  \ 

RHOSL 

Density 

0.02  lbm/in3 

NOTES: 


1  Typical  values  for  4130  steel. 

2  Typical  values  for  FM  5048  silica/phenolic. 

3  Typical  values  for  ATJ  graphite. 

Typical  value  for  phenolic  microballoon. 
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APPENDIX  A 

FORTRAN  LISTING  OF  SIZING  PROGRAM 


A  complete  listing  of  the  FORTRAN  statements  is  furnished  in 
this  appendix  as  a  reference  for  making  any  desired  changes  to  the 
source  deck  of  the  solid  rocket  motor  sizing  program. 
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APPENDIX  B 

DESCRIPTION  OF  MISSILE  OPTIMIZATION  PROGRAM  (MOP* 


B.  1  ANALYSIS 

The  following  analysis*  is  based  primarily  on  the  work  of 
Weisbord  (Ref.  1}  with  modifications  to  simplify  automation  and  no.;  ilo 
changes  to  fit  current  usage. 

The  ideal  impulsive  velocity,  AV];,  achievable  for  any  rockei 
with  constant  specific  impulse  for  each  stage  is  given  in  Equation  1. 

w0l  w02 

AV  =  C!  in  — —  +  C2  in 

WbO! 

where 

WQ.  -  gross  start  weight  for  the  ith  stage,  lbm 
wboi  “  burnout  gross  weight  for  the  ith  stage,  lbm 
Ci  -  gcIsp,  ft/sec 


For  any  given  mission,  AVj  can  be  estimated  from  Equation  2. 


AVj 


*b 


g  sin  0  dt  +  AVd 


(2) 


The  Missile  Optimization  Program  was  developed  by  Mr.  J.  H. 
Dobkins  of  Teledyne  Brown  Engineering. 


where 


AVact  “  actual  impulsive  velocity,  ft/sec 
0  -  angle  of  velocity  vector  above  local  horizon,  rad 

AVj)  -  drag  losses,  ft/sec 

t  -  time,  sec 

t^  -  burn  time,  sec 

For  large  missiles  drag  losses  usually  amount  to  less  than  5  percent 
and  are  not  considered  important  for  preliminary  analyses.  However, 
for  smaller,  high-velocity  vehicles  this  factor  becomes  more  import;:  i. 
and  must  be  estimated. 

The  gravity  loss  term  is  zero  for  horizontal  launch  and  is  (gt)  to. 
a  vertical  launch.  Thus  it  is  strongly  trajectory  dependent  and  some 
approximation  between  the  limits  must  be  made. 

Each  stage  consists  of  motors,  propellant  tanks,  miscellaneous 
hardware,  guidance,  etc.  The  mass  fraction  (XMFj)  and  structure ' 
fraction  (XMSg)  of  stage  i  are  defined  j.n  Equations  3  and  4. 


XMFi  =  Wp./(W0i  -  WPL.)  (3) 

XMSj  =  WS./(W0.  -  WpLi)  (4j 

where 

Wp.  -  weight  of  propellant,  lbm 

Ws^  -  weight  of  stage  -  empty,  lbm 

Wp l.  -  weight  of  stage  payload  (upper  stages 
plus  payload),  lbm. 


Mass  fraction  data  can  be  estimated  from  past  experience.  Structural 
fraction  is  calculated  within  the  program  from  the  mass  fraction  and 
thus  need  not  be  input. 

The  burnout  weight  of  any  stage  may  be  written  in  the  form  of 
Equation  5. 


Wbot  =  W0i+1  +  XMSl  (W0.  -  W0.+1) 


=  XMSiW0i  +  (1  -  XMSi)W0i+1  (5) 


and  Equation  1  becomes 


N  w 

^  „  ”°i 

AVl  =  i?i  C‘  ^  XMS;  W0i  +  (1  -  XMSi)  W0i+1 


(6a) 


or 


n  r  w0.  1 

=  £  Ci  in  [XMS.  Wq.  +  (1  .  XMSO  W0.+J '  AV 


(6b) 


The  problem  requires  that  WD^  be  minimized  for  a  given  AVj. 
Therefore,  the  W0^  become  variables  and  the  maximum  or  minimum 
take-off  weight  is  achieved  when 


3W0l 

8W0.  =  °i=2>N 

°i 


(?) 
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Since  it  is  difficult  to  explicitly  write  these  differentials,  note 
that  Equation  6  is  of  the  form 


£(W0.)  =  0 


(8) 


and  the  partials  can  be  written 


3W0l  3f70WOi 

3W0i  0f/awOl 


(9) 


from  the  chain  rule  for  differentiation.  Equation  9  implies  that  if 
Equation  7  is  to  be  satisfied  then 


9f 


aw, 


=  o  = 


Ol 


3W0; 


(10) 


must  also  be  satisfied.  The  process  yields  for  i  =  2,  3 


3f 


0  = 


C,  {1  -  XMSi  )  C, 


Ci  XMSj 


0WO2  XMS,  w0)  +  {1  -  XMS,  )  W02  W02  XMSi  W0j  +  (1  -  XMSi  )  W0J 


0f 


0  = 


Ci  (1  -  XMSi  )  Cj 


Cj  XMSj 


3W0J  XMSi  W02  +  (1  -  XMSi  )  Wq3  XMSj  w0j  +  (1  -  XMS})  W04 


(11) 


This  can  be  algebraically  manipulated  to  become 


w, 


°3 


C,  XMSi  <1  -  XMS,  ) 


W02  (1  -  XMS,  )  (1  -  XMSj  )  (Cj  -  C,  )  +  Cj  XMS,  (W0j/W02)  (1  -  XMSj  ) 


(12) 


w, 


°4 


Ci  XMS,  (1  -  XMSj  ) 


W0}  (1  -  XMSi  )  (1  -  XMS,  )  (C,  -  Ci  )  +  C,  XMS,  (W02/W03)  (1  -  XMS,  ) 


JB-4 


and  can  be  generally  stated  as 


W0.+1  Ci_i  XMSi  [l  -  XMS  (i  -  1)] 

w0.  (1  -  XMSi.j)  (1  -  XMSj)  (Ci  -  cul)  +  Ci  XMSU1  (Woi  l/W0.)  (1  -  XMSj)  (13) 

i  =  2,N 

Note  that  difficulties  will  occur  when  i  =  1.  An  iterative  process  with 
an  estimated  W02/W0)  ratio  is  required  until  the  required  AVi  is  availa¬ 
ble. 

It  has  not  been  shown  in  the  previous  discussion  that  the  resulting 
data  will  represent  a  minimum  weight  vehicle  (rather  than  maximum). 

A  simple  test  is  to  try  any  other  weight  ratio  between  stages  that  pro¬ 
duces  the  same  performance  and  compare  the  resulting  total  vehicle 
weight,  W0l  .  Other  weight  ratios  invariably  produce  greater  values  of 
W 0z  ,  thus  proving  that  stage  weight  ratios  determined  by  the  previously 
described  procedures  represent  a  minimum  W0l  . 

The  terminal  payload  weight  is  finally  used  to  calculate  the  WQi 
for  the  terminal  stage.  The  vehicle  is  then  specified  totally  from  the 
following  equations. 


Ws«l  ■  Woi+l/(Woi+l/WOi>  -  Woi+l 

(14) 

W_  =  V,'  XMF, 

Pi  Stj  1 

(15) 

wS(  =  wst.  -  Wp. 

(16) 

wboi  =  woi  +  wSi 

(17) 

AVi.  =  C;  in  W0./(W0.  -  Wp.) 

(18) 

where  Wst^  is  the  total  weight  of  stage  i. 


B.  2  PROGRAM  APPLICATION 

The  Missile  Optimization  Program  (MOP)  will  determine  the 
optimum  weight  distribution  between  stages  for  up  to  nine  stages.  If 
for  some  reason  not  enough  stages  were  input,  the  system  will  assume 
extra  stages  with  parameters  equal  to  the  top  stage  until  it  is  possible 
to  build  the  system.  Conversely,  when  the  number  of  stages  is  too 
great  the  last  stage  is  ignored  until  the  system  is  feasible.  This  does 
not  mean  however,  that  the  number  of  stages  has  been  optimized. 

The  input  data  required  is  given  in  Table  B-l.  Provision  has 
been  made  to  study  a  spectrum  of  velocities  and  payloads  and  obtain 
vehicles  optimized  for  each  combination.  Appendix  C  contains  a 
FORTRAN  listing  of  the  program. 


TABLE  B-l .  INPUT  DESCRIPTION  FOR  MOP 


INPUT 

CARD 

TYPE 

VARIABLE 

FORMAT 

1 

DV 

F10.2 

WPL 

FI  0.2 

N 

12 

LAST 

12 

2  (N  cards) 

xisp(r) 

F10.2 

XMF(I) 

FI  0.2 

3 

VST 

F10.2 

WST 

FI  0.2 

JVS 

12 

JWST 

12 

DEFINITION 

Ideal  Velocity,  ft/sec 

Payload,  Ibm 

Number  of  stages,  9  max. 

Last  case  0  No 
1  Yes 

Specific  Impulse,  sec 
Mass  Fraction 

Velocity  Step  Desired 

Weight  Step  Desired 

Number  of  Velocity  Steps  Desired 

Number  of  Weight  Steps  Desired 
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APPENDIX  C 

FORTRAN  LISTING  OF  MISSILE  OPTIMIZATION  PROGRAM  (MOP; 


A  complete  listing  of  the  FORTRAN  statements  is  furnished  in 
this  appendix  as  a  reference  for  making  any  desired  changes  to  the 
source  deck  of  the  missile  optimization  program  (MOP). 


